Conductance quantization in mesoscopic graphene 



N. M. R. Peres, 1 A. H. Castro Neto, 2 and F. Guinea 3 

'Center of Physics and Departamento de Fisica, Universidade do Minho, P-47 10-057, Braga, Portugal 
2 Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA 
3 Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco E28049 Madrid, Spain 

(Dated: February 2, 2008) 

Using a generalized Landauer approach we study the non-linear transport in mesoscopic graphene with zig- 
zag and armchair edges. We find that for clean systems, the low-bias low-temperature conductance, G, of an 
armchair edge system in quantized as G/t = 4ne /h, whereas for a zig-zag edge the quantization changes to 
G/t = 4(n + l/2)e 2 /h, where t is the transmission probability and n is an integer. We also study the effects 
of a non-zero bias, temperature, and magnetic field on the conductance. The magnetic field dependence of the 
quantization plateaus in these systems is somewhat different from the one found in the two-dimensional electron 
gas due to a different Landau level quantization. 

PACS numbers: 71.10.Fd, 73.23.Ad, 72.10.Bg, 72.I5.Lh. 



I. INTRODUCTION 

Graphene, a two-dimensional (2D) carbon system on a 
honeycomb lattice presents many anomalous properties when 
compared with the well-known 2D electron gas obtained 
in heterostructures. One of the most striking properties 
is an exotic integer quantum Hall effect (IQHE) predicted 
theoretically 1,2 , and measured recently 3 - 4 . The IQHE shows 
a Hall conductivity given by: a xy = 2(2n + l)e 2 /h, where 
n is a positive integer. Interestingly, the electrical properties 
of the graphene systems can be considered classical, in the 
sense that the measured conductance of the systems is found 
to increase with the increase of system width and to decrease 
with the increase of the system length 5 . This experimental re- 
sult can be understood as an evidence for the presence of dis- 
order in the measured systems. This is further supported by 
the difficulty in finding experimental evidence for a fractional 
quantum Hall effect (FQHE)i^. 

In the closely related field of carbon nanotubes, recent ex- 
periments showed that the conductance of a single wall car- 
bon nanotube is quantized 7 and shows Fabry-Perot interfer- 
ence patterns. These results can be explained within a gen- 
eralized Landauer approach, S— matrix theory^, and non- 
equilibrium Green's function methods 910 . The formulation 
of the problem was introduced by Lake et al. 1 1 , after the work 
of Caroli et Because carbon nanotubes are essen- 

tially wrapped graphene, we expect conductance quantization 
and Fabry-Perot interference patterns to be also observable in 
ultra-clean graphene. The quantization and the interference 
patterns, however, should reflect the different types of edges a 
graphene sheet has (see Fig.^. 

The importance of zig-zag and armchair edges in graphene 
sheets has been recognized in electron microscopy 15 . Clearly, 
these two types of edges produce very different electron mi- 
croscopy intensity curves. We expect, therefore, that coherent 
charge transport should be different if measured in systems 
with different types of edges. As in carbon nanotubes^, a 
simple Landauer approach to determine the quantization of 
the conductance, G, can be used for mesoscopic graphene 
sheets. The calculations follow the generalization of the Lan- 



dauer approach introduced earlier by Bagwell and Orlando 17 . 
This type of approach does not account for a discussion of 
interference patterns, since it neglects multiple electronic re- 
flection at the contacts 7,9 , and it will be discussed elsewhere 18 . 
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FIG. 1: (color online) Geometry of a finite-size honeycomb lattice. 
Top: Sample with a zig-zag edge; Bottom: Sample with an armchair 
edge. The thick dotted lines represent the position of the leads which 
are assumed to be made out of graphene. The rectangles in the top 
left corner of the figures (close to the left lead) show the geometry of 
the unit cell. 

This paper is organized as follows: in Sec. [H] we introduce 
the tight-binding solution for electrons in graphene strips hav- 
ing two different types of edges and the differences and simi- 
larities between graphene and carbon nanotubes are discussed. 
Starting from the tight binding solution of laterally confined 
electron in graphene strips we introduce a Landauer analysis 
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allowing for the calculation of the conductance due to quasi- 
one-dimensional transport arising from the confinement; this 
is done at at zero-temperature and in the zero-bias limit. In 
Sec. [HI] we discuss the effect of temperature, gate bias, and 
magnetic field on the conductance curves. It is found that con- 
ductance has a minimum value as function of the gate poten- 
tial (or Fermi energy), from which a "V-like" curve develops, 
in agreement with measured transport properties in graphene 
samples. Sec. lIVI contains our conclusions. 



IL A LANDAUER ANALYSIS 

The geometry of the problem is shown in Fig.Q where the 
left and right leads are supposed to be made of graphene 9 . 
This aspect is not essential in the approach we develop below, 
where the contacts to the leads will be modeled by a trans- 
mission probability (as done in Refs. II7I9I1 in the context of 
carbon nanotubes). The systems studied here are assumed to 
have a very asymmetric aspect ratio, where the length, L, is 
much larger than their width, W. These systems have some 
similarity with carbon nanotubes 16 but differ from them in a 
fundamental way: the absence of periodic boundary condi- 
tions along the direction perpendicular to the edges (the m 
direction). As a consequence, it is possible to have different 
kinds of strips, characterized by different types of edges. In 
what follows, we discuss the cases of zig-zag and armchair 
edges although other edge geometries can be studied with the 
same methods. 

The calculation of the conductance G following a Landauer 
type of approach 1617 requires the solution of a tight-binding 
problem in a finite geometry. The tight-binding Hamiltonian 
has the form : 

Ht.h. = -t J2 (4 a* + kc -) 

+ t' J2 ( a U%>+ & l,A- +h - c -)' (!) 

where a\ a (ai jCr ) creates (annihilates) an electron on site R;, 

with spin a (a =|, J,) on sub-lattice A and b\ a (bi.a) creates 
(annihilates) and electron on site Rj with spin a (a =\ , J.) on 
sub-lattice B. t is the nearest neighbor ((i, j)) hopping energy 
(i w 2.7 eV), and t' is the next-nearest neighbor (((«, j))) 
hopping energy (f /t « 0.1). In what follows we suppress 
the spin index since it plays no role (apart from a degeneracy 
factor). 

In an infinite system the Hamiltonian can be easily di- 
agonalized and one can show that the low energy electronic 
excitations of the problem reside around the K-points of the 
Brillouin zone 19 and have a dispersion given by (t' does not 
remove the Dirac point): 



E±(k) =±v F \k\ 



(2) 



where k = (k x ,k y ) is a two-dimensional momentum, and 
vp = 3ta/(2h) (where a is the lattice spacing) is the Fermi- 
Dirac velocity. Eq. (|2ji describes the dispersion relation of 



Dirac electrons with a speed v p ■ One of the consequences of 
the Dirac dispersion is that the fermions in the system have 
zero effective mass, and a linearly vanishing density of states, 
N(E) (N(E) oc \E\), at low energies. The linearly dispers- 
ing electrons and the vanishing of the density of states lead 
to a very anomalous metallic behavior with many non-Fermi 
liquid properties 1 . These anomalous properties are reflected 
in the experimentally measurable quantities, such as the Hall 
conductivity in the IQHE^i. We are going to show that the 
presence of Dirac fermions in the spectrum has also a strong 
effect in the conductance of finite graphene strips. 
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FIG. 2: (color online) Label of the carbon atoms for Hamiltonian 
Q in the bulk of the system (only the case t' = is shown for 
simplicity). The reader is referred to Fig. Q where the unit cell is 
shown in each case. The black circles, labeled a, b, c, and d, inside 
the rectangle or on the rectangle edges, belong to same unit cell; the 
other circles represent lattice points in adjacent unit cells connect by 
the hopping matrix t. 
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FIG. 3: (color online) This figure shows how to deal with the tight- 
binding boundary conditions introduced by the free edges of the sam- 
ple. The cells on the right refer to top free edge and the bottom free 
edge of a zig-zag sample (top of Fig.Q). The cells on the left refer 
to the top free edge and the bottom free edge of an arm-chair sample 
(bottom of Fig.Q). As in Fig.|2| only the sites connected by t are 
shown. 
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In a finite system the boundary conditions (and hence, the 
type of edges) become important in order to define the energy 
spectrum. In Fig.|2]we establish the labeling of the carbon 
sites in the bulk of the honeycomb lattice, and in Fig. [5] we 
present the labels for the carbon atoms close to the free edges 
of the system. Notice that each unit cell in Fig. |3 and Fig. [3] 
contains four carbon atoms, labeled by wave function ampli- 
tudes a(n, m), b(n, m), c(n, m), and d(n, m), where n and m 
are integers that label each unit cell. Since we are assuming 
that in each system only the edges are different, we impose 
periodic boundary conditions along the direction parallel to 
the edges (the n direction,) leading to one-dimensional (ID) 
transport where the electronic states can be labeled by the mo- 
mentum q x . If N y is the number of unit cells along the m di- 
rection, then the the tight-binding problem with zig-zag edges 
has a dimension given by (4N y + 2) x (AN y + 2), while for the 
armchair edge its dimension is (4^ + 4) x (AN y + 4). For 
a nanotube, with periodic boundary conditions, the problem 
can be formulated in terms of two amplitudes instead of four. 

The calculation of the conductance of a two-dimensional 
lattice system with a very asymmetric aspect ratic 16 requires 
the identification of the number of transverse modes, M(e), 
at a given energy, e. M (e) can be obtained from the solu- 
tion of the corresponding tight-binding problem. In Fig. 0] we 
show the energy bands obtained from the diagonalization of 
the tight-binding Hamiltonian Q with t' = for the zig-zag 
and armchair edges. It is clear that the two different types 
of edges lead to two very different band structures, especially 
close to zero energy. 
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FIG. 4: (color online) Energy levels for t' = (in units of t) as 
function of q x (in units of the lattice spacing along the n direction) 
for the two different types of edges. N y is the number of unit cells 
perpendicular to the edges (along the m direction) and N x gives the 
number of momenta used in the abscissa (it also corresponds to the 
number of unit cells along the n direction). The horizontal line at 
zero energy on the right panels has an interception with zero energy 
modes in the upper panel (metallic behavior) and no interception in 
the lower panel (semi-conducting behavior). 



It is important to mention the similarities and the differ- 
ences between planar systems 20 and carbon nanotubes— . An 
armchair nanotube has the hexagons having two sides perpen- 
dicular to the tube axis, in the zig-zag nanotube the hexagons 
have two sides parallel to the tube axis. A graphene sam- 
ple with a zig-zag edge has an energy spectrum presenting 
some similarities with an armchair nanotube. It has two bands 
crossing the chemical potential at zero energy (and finite mo- 
mentum). A graphene sample with an armchair edge has a gap 
at zero energy (near q x = 0) and hence is insulating, except 
when N y , the number of unit cells in the m-direction, is a mul- 
tiple of 3, in which case the gap goes to zero and the material 
is metallic (in fact a zero gap semi-conductor; see Fig.[4j 20 . 

A generalized Landauer approach 17 shows that the tunnel- 
ing current is given by: 



I(V,T) 



2c 
17 



d e M(e)t(e,y)[/(e-/ii)-/(e-/i 2 )], (3) 



where t(e, V) is the transmission probability per conducting 
mode at the energy e, /(e) = 1/ (e e/T + 1) is the Fermi-Dirac 
distribution (T is the temperature and we have put ks = 1), 
V is the bias voltage applied to the system, and /ii (112) is the 
chemical potential at right (left) lead {jii =112 + eV). 

In a clean system, all ID modes can carry electric current, 
as long as they have a finite velocity in the direction parallel 
to the edge. Although the zig-zag edge system has zero en- 
ergy modes with finite q x momentum, the group velocity of 
these modes is zero when t' = 0, and therefore they do not 
contribute to the conductance. If we neglect the effect of the 
next nearest neighbor hopping (t' — 0), both edge systems 
have two conducting zero energy modes (choosing a metallic 
armchair edge system-iY y multiple of three), and as a conse- 
quence the small bias conductance is given by 4e 2 t/h, at zero 
chemical potential. 

With the addition of a next nearest neighbor hopping (t' 7^ 
0) the picture for the zig-zag edge changes substantially: (1) 
the half-filling case occurs at finite energy (not at zero energy); 
(2) the flat band states located at zero energy, for t' = 0, be- 
come dispersive for t' 7^ 0, and therefore the conductance is 
modified. Fig. |5] shows the energy levels as function of the 
momentum q x for a zig-zag system with t' = and t' = Q.2t, 
for N y = 3 (this rather small N y value allows the individual 
visualization of the transverse modes over the full bandwidth, 
as it does in Fig.|4]as well). 

In order to determine the conductance it is necessary to 
know how many ID transverse modes are active for transport 
at a given energy. For an armchair edge with N y multiple of 3 
(and hence a metallic system) one has two conducting modes 
at zero energy (see Fig. |6} . As one moves away from zero 
energy the number of modes increases as shown in Fig. |S] 
Hence, the zero-bias zero-temperature conductance of a clean 
metallic armchair system is given by: 



(4) 



with n a positive integer. The value of n depends on the value 
of the gate potential controlling the electronic density. The 
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FIG. 5: (color online) Energy levels (in units of t) as function of 
q x (in units of the lattice spacing along the n (or x) direction) for a 
zig-zag system. Right: t' = 0; left: t' = Q.2t 



conductance also depends on t, the transmission probability, 
assumed to be energy independent for simplicity. 

A zig-zag clean system with t' = shows a different de- 
pendence. At zero energy G is given by Ae 2 t/h because of 
the presence of the two conducting zero-energy modes (see 
Fig. [6}. However, as the gate potential moves slightly away 
from zero the value of the conductance drops to 2e 2 i/h, since 
only one transverse mode is available. So the situation with 
a zig-zag system is somewhat ill-defined. It reasonable to ex- 
pect that this sudden change on the conductance is difficult 
to be experimentally observed, since it would require an ex- 
treme fine tuning of the experimental parameters. Therefore, 
we expect that the zero-bias zero-temperature conductance of 
a system with a zig-zag edge to be given by: 



G 



2(2n+ l)e 2 ~ 



t, 



(5) 



with n depending also on the value of the gate potential. 

The special value of G we found exactly at zero energy (in 
the zig-zag case) does not survive, however, when electron- 
hole symmetry (by this it is understood that the energy spec- 
trum is not symmetric around zero energy) is broken by the 
presence of a finite t'. When t' ^ the zero mode acquires 
dispersion and the conductance is given by Eq. (0. In Fig. [7] 
the effect of t' in the low-energy band structure and in M(e) 
is visible. Clearly the zero energy modes have been removed, 
leading to two degenerate dispersive bands, and an asymme- 
try in energy in the steps of M is introduced. The asymmetry 
in energy of M(e) is an experimental way of measuring the 
value of t' in these systems. In the armchair case, the effect 
of t' on G is not as dramatic as in the case of the zig-zag edge 
system, since in this case the conductance does not have an 
abrupt change around zero energy. 

In a system with several graphene planes each plane con- 
tributes to the conductance almost independently because of 



the weak coupling between graphene sheets. Therefore, the 
resulting conductance should be, at least approximately, given 
by the above results multiplied by the number of layers. 
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FIG. 6: (color online) Lower panels: one-dimensional energy bands 
(energy in units of t) associated with a zig-zag (left) and a armchair 
edge (right) systems for t' = 0. Upper panels: number of ID chan- 
nels, M, as a function of energy (in units of t). 
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FIG. 7: (color online) Left panel: one-dimensional energy band (en- 
ergy in units of t), for a zig-zag edge system with t' = 0.2t. Right 
panel: number of ID channels, M, as a function of energy (in units 
of t). The horizontal dashed line helps in stressing the fact that M(e) 
is asymmetric as a consequence of a finite t'. 



III. TEMPERATURE, MAGNETIC FIELD, AND VOLTAGE 
BIAS EFFECTS. 

In this section we discuss how the zero-bias and zero- 
temperature results are modified by considering the more gen- 
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eral case of a finite bias, temperature, and magnetic field ap- 
plied perpendicular to the graphene plane. 

Close to equilibrium (where Hi — H2 = Ep) the conduc- 
tance can be obtained from (|3} as: 



G(V 



2e 2 



deM(e)t(E, V) 



df(e - Ep) 
de 



(6) 

For finite bias the conductance is determined from Eq. (0 
after a simple numerical derivative in relative to the bias po- 
tential. Notice that in equilibrium, the changes in Ep can 
be obtained by simple changes in the value of the gate volt- 
age. The results for G(V, T) as function of Ep, V, and T, are 
shown in Fig. [8] One can clearly sees that, as predicted, the 
conductance is quantized in units of 2te 2 /h, being even in the 
case of the armchair edge and odd in the case of the zig-zag 
edge. The temperature makes the plateaus in the conductance 
smooth. Application of a bias voltage V shifts the position of 
the conductance plateaus, as expected. 
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FIG. 8: (color online) Conductance (in units of 2te 2 /h) for armchair 
and zig-zag systems as function of Ef, for various values of T and 
V(t = 2.7eV). 

The effect of an external magnetic field in the conductance 
quantization of the 2D electron gas was experimentally stud- 
ied in the past 23 and discussed in general terms by Btittiker— . 
We consider only the case of systems with zig-zag edges for 
simplicity since, with periodic boundary conditions along the 
n direction, a unit cell with only two atoms can be chosen. In 
the presence of an applied perpendicular magnetic field B the 
hopping integrals change to: 



U 



tijC 



(7) 



where the phases 0jj = 27r J? A ■ dl/tfro, and </>o = h/e is 
the quantum of flux. The sum of 9ij over a close path defined 
by the hopping integrals must equal the value 27r0p/0o = 
2irBAp/(fir), where Ap is the area enclosed by the path P. 
For the honeycomb lattice, the equation for the amplitudes of 



the tight-binding Hamiltonian Q can be written as: 



— t\h 4- h -L e i27r <i >m h 1 

^[yn^m T Un— 1,771 1 c u n,m— 1J 

,/r i7T0/3 , -i27T<A(lTi-5/6) 

-427T0(m-l/6) , -i7T0/3 

i2^(m+l-5/6) , i27T0(m-l/6) 1 



(8) 



+ 6„-i, m +ie l2 ^ (m+1 - 1/6) +e a ^ m - 5 ^b n ^ m ] , 

(9) 

where <\> = _B^4 c /</)o, and A c is the area of an hexagon (A c — 
3v / 3a 2 /2 with aw 1.4 A in graphene). For <fi — we obtain 
the results of Sec.llTI 

In the presence of a magnetic field the states of the bulk 
graphene are described in terms of Landau levels. At low en- 
ergies, when the Dirac fermion description © is valid, the 
energy levels are given by: 



E±(n) 



2a 



2 7 2 



(10) 



where we have assumed t' <§; t, and defined lp — \J h/eB 
as the magnetic length, a — 9t'a 2 /4, and 7 = 3ta/2 
in = 1, 2, . . . ). For t' = 0, the energy levels are given by: 
E±(n) — ±\/2^l E ^ \fn. This result shows that, for the case 
of Dirac fermions, and unlike the ordinary 2D electron gas, 
the Landau levels are not equally spaced 2 ^. Notice that the 
cyclotron energy, fku c = ^/2vph/lB, is much larger than the 
Zeeman energy, gfipB ig sa 2, fip is the Bohr magneton - 
for B = 12 T, hw c fa 0.142 eV and g/j, B B w 7 x 10" 4 eV). 
Thus, we disregard the Zeeman energy in what follows. 

In a finite system the energy levels given by i ll Ob are modi- 
fied by the lattice structure and by the presence of edges. This 
can be clearly seen in Fig. [9] where we plot the solution of (|9jl 
for a graphene strip with a zig-zag edge. Clearly at B = 10 T 
the Landau level structure predicted by (II Oi is seen to develop 
close to the Dirac point. The effect of the magnetic field on 
G is two fold: (1) the magnetic field leads to non-dispersive 
magnetic levels, with a large degeneracy; (2) the energy level 
spacing is modified giving rise to a piling up of energy levels 
as one moves away from the Dirac point. 

Notice that our discussion is valid for weak magnetic fields 
and hence does not apply to the quantum Hall regime that was 
discussed in refs. II lO . We observe that states with values 
of q x away from the system edges can not contribute to the 
conductance, since their group velocity is zero. Only those 
levels having non-zero group velocity can act as ID channels 
for electron transport. The piling up of energy levels has the 
experimental consequence that the observation of many quan- 
tized plateaus becomes difficult. In addition, the increase of 
the degeneracy of each Landau level with the increase of the 
field reduces the number of observable plateaus, as in the nor- 
mal electron gas 23 . On the one hand, when we compare M(e) 
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FIG. 9: (color online) Energy spectrum (energy in units of t) of a 
graphene strip with a zig-zag edge in a magnetic field B (Left: B = 
10 T; Right: B = IT), Lower left panel: number ID transverse 
modes M for B = 10 T. 



for B = and i? = 10 T we see that energy width of the 
AI = 1 step in the latter case has a much larger value. On 
the other hand, the piling up of the Landau levels leads to the 
reduction the energy value of plateaus (see Figs. Qand|9}. A 
small field B does not lead the formation of Landau levels, but 
removes the degeneracy of the band formed from the flat band 
of zero modes when t' is considered, and leads to a M = 1 
energy step width of a larger value when compared with the 
B = case. 



point, and therefore the conductance assumes a "V"-shape 
form as a function of the gate potential. The lowest value 
of G/t in the zig-zag edge system is 2e 2 /h (for a non-zero, 
albeit small, t'), whereas in an armchair edge system we find 
4e 2 /h. 

We have studied in detail the plateaus in the conductance 
of graphene strips as a function of temperature, applied gate 
voltage, and external magnetic field. We have seen that the 
temperature smoothes out the plateaus and that applied gate 
voltages shift the plateaus in energy. We have also discussed 
the effect of next nearest neighbor hopping energy t', that 
breaks the particle-hole symmetry of the problem and intro- 
duces dispersion for in the zero modes. The effect of a finite 
magnetic field is quite interesting in these systems because of 
the unusual relation between the energy and the Landau level 
index. We show that a magnetic field has effect in piling up 
the conductance steps and modify their size in energy. These 
effects should be easily observable in ultra clean mesoscopic 
graphene strips. 

For graphene samples of 10-100 fim size 25 it was found 
that the conductivity, given by a = GL/W, has the univer- 
sal value of ~ 4e 2 //i. This result can be understood using 
a bulk calculation of the effect of vacancies on the electric 
linear response (Kubo formula) of Dirac fermions 1 . ^From 
the point of view of coherent tunneling, these experimental 
results indicate that these samples are in the ohmic regime, 
having a mean free path shorter than the system size. We be- 
lieve, however, that in ultra clean graphene samples it will be 
possible to observe conductance quantization and interference 
patterns, as it is the case of carbon nanotubes. We hope that 
our results will stimulate further studies of transport in these 
amazing systems. 



IV. SUMMARY AND CONCLUDING REMARKS 

We have discussed the tunneling transport in clean meso- 
scopic graphene strips. We show that different graphene strips 
have different conductance values due to different types of 
edges. As a general consequence of the graphene band struc- 
ture, and at odds to the usual 2D electron gas systems, the 
conductance always increase as we move away from the Dirac 



Acknowledgments 

We thank A. Geim, P. Kim, and W. de Heer for stim- 
ulating discussions. N.M.R.P. thanks ESF Science Pro- 
gramme INSTANS 2005-2010 and FCT under the grant 
POCTI/FIS/58 133/2004. A.H.C.N. was supported through 
NSF grant DMR-0343790. 



1 N. M. R. Peres, F. Guinea, and A. H. Castro Neto 
cond-mat/0506709 and cond-mat/05 12091 

2 V. P. Gusy nin, and S. G. Sharapov, Phys. Rev. Lett. 95, 146801 
(2005), and cond-mat/05 12 157 

3 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. Kat- 
snelson, I.V. Grigorieva, S.V. Dubonos, and A.A. Firsov, Nature 
438, 197 (2005). 

4 Yuanbo Zhang, Yan-Wen Tan, Horst L. Stormer, and Philip Kim, 
Nature 438, 201 (2005). 

Claire Berger, Zhimin Song, Tianbo Li, Xuebin Li, Asmerom Y. 
Ogbazghi, Rui Feng, Zhenting Dai, Alexei N. Marchenkov, Ed- 
ward H. Conrad, Phillip N. First, and Walt A. de Heer, J. Phys. 
Chem., 108, 19912 (2004). 



6 A. H. Castro Neto, F. Guinea, and N. M. R. Peres, 
cond-mat/0509709 

7 Wenjie Liang, Marc Bockrath, Dolores Bozovic, Jason H. Hafner, 
M. Tinkham, and Hongkun Park, Nature 411, 665 (2001). 

8 For a compreensive presentation of the Landauer formulation 
see Supriyo Datta, Electronic Transport in Mesoscopic Systems, 
(Cambridge, 1997). 

9 S. Krompiewski, J. Martinek, and J. Barnas, Phys. Rev. B 66, 
073412 (2002). 

10 For a review on the non-quilibrium Green's function method see: 
J. Rammer, and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

11 Roger Lake, Gerhard Klimeck, R. Chris Bowen, and Dejan Jo- 
vanovic, J. Appl. Phys. 81,7845 (1997). 



7 



C. Caroli, R. Combescot, P. Nozieres, and D. Saint- James, J. Phys. 
CSolid St. Phys., 4, 916 (1971); idem, ibidem, 4, 2598 (1971); 
idem, ibidem, 5, 21 (1972). 

13 R. Combescot, J. Phys. CSolid St. Phys., 4, 2611 (1971). 

14 R. Combescot, and G. Schreder, J. Phys. C:Solid St. Phys., 6, 1363 
(1973); idem, ibidem, 1, 1318 (1974). 

15 Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. Tsukada, and 
Hiroshi Fukuyama, Appl. Surf. Sci. 241, 43 (2005). 

16 Weidong Tian, and Supriyo Datta, Phys. Rev. B 49, 5097 (1994). 

17 Philip F. Bagwell, and Terry P. Orlando, Phys. Rev. B 40, 1456 
(1989). 

18 N. M. R. Peres, F. Guinea, and A. H. Castro Neto, in preparation. 

19 P. R. Wallace, Phys. Rev. 47, 622 (1947). 



Kyoko Nakada, Mitsutaka Fujita, Gene Dresselhaus, and Mildred 
S. Dresselhaus, Phys. Rev. B 54, 17954 (1996). 
Riichiro Saito, Mitsutaka Fujita, G. Dresselhaus, and M. S. Dres- 
selhaus, Phy. Rev. B 46, 1804 (1992). 

M. H. Johnson and B. A. Lippmann, Phys. Rev. 76, 828 (1949); I. 
I. Rabi, Z. Phys. 49, 507 (1928). 

B. J. van Wees, L. P. Kouwenhoven, H. van Houten, C. W. J. 
Beenakker, J. E. Mooij, C. T. Foxon, and J. J. Harris, Phys. Rev. 
B 38, 3625 (1988). 

M. Buttiker, Phys. Rev. B 38, 9375 (1988). 

K. S. Novoselov, D. Jiang, T. Booth, V.V. Khotkevich, S. M. Mo- 

rozov, and A. K. Geim, Proc. Natl. Acad. Sci. 102, 10451 (2005). 



